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Abstract. Using quantum theory operator methods we discuss the general reversible reactions A\ + A2 + ■ ■ ■ A T <-> 
Bi + B2 + ■ ■ ■ + B s , where r and s are arbitrary natural positive numbers. We show that if either direction of the 
reaction is repeated a large number of times N in a finite total time T then in the limit of very large N , keeping 
T constant, one remains with the initial reacting particles only. We also show that if the reaction evolves through 
different possible paths of evolution, each of them beginning at the same side of the reaction, proceeds through different 
intermediate consecutive reactions and ends at the other side, then one may "realize" any such path by performing 
in a dense manner the set of reactions along it. The same results are also numerically demonstrated for the specific 
reversible reaction A + B <-> A + C. We note that similar results have been shown to hold also in the quantum regime. 
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1. Introduction 



The general reversible reactions A\ + A2 + • • • + A r <-> B\ + B2 + • • • + B s , where r and s are two 
arbitrary positive natural numbers, have been studied by many authors (see, for example, (Ben- 
Avraham and Havlin, 2000) and references therein). These studies discuss, especially, the effects of 
the single reaction, or, in case it is repeated N times, the effect of these repetitions where the general 
total time T increases proportionally to N. We can, however, imagine a situation in which the rate 
of these repetitions increases and discuss the effect of this increase upon the reaction. Such an effect 
have been studied in (Kac, 1947) with respect to random walk and it was shown that when the rate 
of repeating this walk becomes very large one obtains a Brownian motion (Nelson, 1967). It has also 
been shown (Bar, 2001), (Bar, 2003) with respect to the finite one-dimensional array of N imperfect 
traps (Ben-Avraham and Havlin, 2000), (Smoluchowski, 1917) that as N becomes very large the 
survival probability (Ben-Avraham and Havlin, 2000) of the particles that pass through them tends 
to unity. We show in this work that if either direction of the mentioned reversible reaction is repeated 
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a large number N of times in a finite total time T, then in the limit of very large N (where T is 
kept constant) one remains with the initial reacting particles only. 

We use quantum theory methods and terminology as done by many authors that use quantum 
formalism for analysing classical reactions (see for example (Mattis and Glasser, 1998) and annotated 
bibligraphy therein). The most suitable quantum method is the coherent state one (Glauber, 1963), 
(Klauder and Sudarshan, 1968) which allows us to define simultaneously, as has been remarked in 
(Swanson, 1992), the expectation values of the conjugate variables Q and P. Thus, they may both 
have nonzero values in which case the entailed formalism resembles (Swanson, 1992) the classical 
one and so it may be used for discussing classical reactions. 

The use of the coherent state formalism, together with second quantization methods, for classical 
systems have been extensively studied by Masao (Masao, 1976). Since the described phenomena and, 
especially, the particles participating in the reactions are classical we represent them here by real 
coherent states. That is, we denote the reacting and product particles by the real coherent states 
\z > = e 2 1^1 J2nZn^ z i \ n >, where z is the real number z = q and p are two arbitrary 

(n!)2 (2)3 

real c numbers and \n > are number representation eigenstates (Klauder and Sudarshan, 1968). 
The masses of the reacting and product particles are assigned, for convenience, the unity value. 
We note that although the mathematical entities and "operators" involved in this method do not 
conform, as will be shown, to the known quantum operator formalism, we, nevertheless, follow, 
except for the differences, the conventional definitions and methods of the last theory. The results 
obtained by applying the real coherent state formalism for classical reactions are exactly the same as 
those previously obtained (Bar, 2000) by applying the complex coherent state methods for quantum 
reactions. Moreover, although the real coherent state formalism entails an evolution operator (see 
the following Eq (3)) which is nonunitary and unbounded we, nevertheless, show that the results 
obtained are valid also for this kind of operator. That is, we obtain for the classical reactions the 
same results that were obtained (Bar, 2000) for the quantum analogues which are discussed by using 
complex coherent state formalism (Glauber, 1963), (Klauder and Sudarshan, 1968) which involves 
unitary and bounded evolution operators. Needless to say that one may use the conventional complex 
coherent state formalism (Glauber, 1963), (Klauder and Sudarshan, 1968) for discussing classical 
reactions. 

In Section 2 we discuss the reversible reaction A\ + A2 + • • • + A r <-> B\ + B2 + ■ ■ ■ + B s for the 
special cases of r = s = 1 and r = s = 2. That is, we study the reversible reactions A <-> B and 
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A-y + A 2 <-> Bi + B 2 , and show that repeating either direction of each a large number N of times in 
a finite total time T results, in the limit of very large N, in a unity probability to remain with the 
initial reacting particles only. The generalization to any finite r and s follows. 

We note that since we discuss the probability to remain with the initial reacting particles the 
product particles of such reactions are not relevant (as the reacting ones) to our discussion. In 
Section 3 we discuss the more general and natural case in which the product particles are relevant. 
That is, we assume that the particles of the ensemble interact at different places and times and 
that they begin from some given initial configuration of reactions and end at a final different one. 
We assume that there are large number of different paths of reactions that all begin at the given 
initial configuration and end at the final one and we calculate the probability that a specified 
system of reacting particles evolves along some prescribed path of them. We note that such paths 
of "states" for the diffusion controlled reactions have been discussed in (Masao, 1976), (Mikhailov, 
1981), (Mikhailov and Yashin, 1985), (Namiki, 1992) where use was made of quantum field theory 
methods (Mahan, 1993), (Enz, 1992), (Mattuck, 1976), including Wick's theorem (Mahan, 1993), 
(Enz, 1992), (Mattuck, 1976), to derive the classical Feynman diagrams and paths (Feynman, 1948), 
(Feynman and Hibbs, 1965). These methods were also used in (Mikhailov, 1981), (Mikhailov and 
Yashin, 1985) for chemical kinetics. We show that taking the limit of a very large number N of 
reactions along the prescribed path and performing them in a dense manner one obtains a unity 
value for the probability of evolution along that path. 

In Section 4 we use a numerical model that has been used in (Bar, 2001) for showing the effect 
of dense measurement on classical systems. This is the model of two dimensional concentric billiard 
(Bar, 2001) that is used here to numerically simulate the reversible reaction A + B <-> A + C. 
The two possible modes of reflections inside the billiard; either between the two concentric circles 
or between points of the outer circle, represent the two directions of the reaction. We note that 
nuclear and radioactive reactions are well simulated by billiards in which the stationary scattering 
circles represent the interactions between particles (see, for example, (Bauer and Bertsch, 1990) in 
which a model of a rectangular billiard with a circle inside was used to discuss the decay law of 
classical systems). We show that the numerical billiard simulations conform to the analytical results 
of Sections 2 and 3. 
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We assume that we have some set of N identical particles so that the configuration in which the 
ith particle is located at qi (i = 1,2, ■■ ■ , N) is defined as a state of the system and denoted, in 
Dirac's notation, \q±,q2, • • • ,qN > {\q N >) (for distinguished sets of particles we partition (Masao, 
1976) the total N system to N±, N 2 , ■■■ subsystems). Thus, when representing, in the following, 
classical particles by states we mean that they are elements of some configuration of the whole 
system. Following this terminology we may calculate the probability to find the set of particles in 
some definite state \q N > as (Masao, 1976) 

F (N Hqi,q2, ■■■,q N ;t)= £ f< N \ Ql , q 2 , ■■■,q N ;t), 

all permutations of qi 

where f^ N \qi, q2, • • • , QN', t) is some normalized distribution function. To this probability one as- 
sign, as done in (Masao, 1976), a "state" \F(t)>= J2n=o I dQ N F^ N \q±,q 2 , ■ ■ ■ ,qN',t)\q N >, where 
SiQ N = LgL (the divlslon by m K necessaIy (M asao, 1976) so as not to overcount the state 
\q N > AM times). Thus, the former probability to find the system in the state \q N > may be written 
as (Masao, 1976) F (N \ qi ,q 2 , ■ ■ ■ ,q N ;t) =<q N \F( y t)>. 

We discuss first the general reversible reaction A\ + A 2 + • • • A r <-> B\ + B 2 + • • • + B s for the 
specific case of r = s = 1, A <^ B where A and B are, as noted, represented by the two coherent 
states (Klauder and Sudarshan, 1968) 

n=oo n 

\z A >=e~W £ 3*|n> (1) 

n=0 (™0 2 
n=oo n 

\ ZB>=e -\\*B? £ ^|n> 

n=0 (™0 2 

Using the following general equation for any two operators X and Y (Klauder and Sudarshan, 1968) 

e Y Xe~ Y = X + [Y,X} + ±[Y, [Y, X}} + • • • , 
where [Y, X] is the commutation [Y, X] = YX — XY, one obtains 



U(q,p)(aP + PQ)U-\q,p) = a(P + p) + (3{Q + q) 



(2) 
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The a, (5 are arbitrary parameters, U(q,p) and U~ 1 (q,p) are given respectively by U(q,p) = e p ®~ qP , 
U^ 1 (q,p) = U(—q,—p), and Q, P are the coherent state operators that satisfy [Qi,Pj] = 5ij. That 
is, U(q,p) translates the operators Q and P by q and p respectively. Now, since the coherent states 
have been defined, as remarked, in terms of the number representation eigenstates (see Eq (1)) we 
write the time evolution operator of the relevant states as e Nt , where N is the number operator 
(Klauder and Sudarshan, 1968) (note that since we discuss in this paper real coherent states the 
evolution operator is real also). 



N is defined analogously to the corresponding operator of the complex coherent state formalism 
(Glauber, 1963), (Klauder and Sudarshan, 1968) but without the complex notation i in the middle 
expression. Note that the operator N is not positive definite and this to remind us, as remarked, 
that the real coherent state formalism discussed here does not conform to the conventional quantum 
operator process. Nevertheless, as remarked, the final results obtained here are exactly the same as 
those accepted in (Bar, 2000) from the quantum complex coherent state formalism. The commuta- 
tion [Qi,Pj] = Sij have been used at the right hand side of the last equation. Applying the operator 
TV on the coherent state \z > from Eq (1), and taking the scalar product of the result with the 
conjugate state < z\ one obtain (using < n\e Nt \m >= e nt 5 nm , since in the number representation 
the operator N is diagonal) 



N = a) a = (9^ 7 Z ){ 9±Z ) = I ( q2 _ p 2 1} 



<z\e \z>= exp( 




^l 2 )E 



ze l z 



) 



,n 



n=0 



n! 



= exp(— -|z| 2 — -^\z\ 2 + ze t z) =<z\e t z> 
=<z|(cosht + smh.t)z>=<q,p\q t ,pt > 



(3) 



The last result is obtained by writing z in terms of q,p in which we have 



qt = g(cosht + sinht), p t = p(cosht + sinht) 



(4) 



We, now, calculate, using Eq (3), the probability p(\qAiPA>) to remain with the initial particle A 
after the reaction A — > B where the particle B is represented by the coherent state e Nt \zA>- This 



The effect of increasing the rate. 



7 



is given by 




iil^il — \j — \ • • " -y 

= exp(--(p A + ?a) - ^(pa, + q At ) ) 2^ 

n=0 



= exp(--(pA + 9a) 2 - ^(pA t + qA t f + +PA) • 

■{qA t +PA t )) = exp(-^(pA + QA) 2 - ^(PA + 9a ) 2 (cosh t + sinh t) 2 + 

+ ^(qa + Pa) 2 ■ (cosht + sinht)) = exp(--{p A + Qa) 2 (^ + 

+ ^(cosht + sinhi) 2 — (cosht + sinht))) 



(5) 



Note that since we discuss coherent states the interpretation (Klauder and Sudarshan, 1968) of the 
expression < ZA\e Nt \za > is the probability to find the mean position and momentum of the coherent 
state e Nt \zA>, which represents B, equal to those of za, which represents A, and this probability 
is equivalent in our discussion to remaining with the particle A. Prom Eq (5) one obtains the 
probability to remain with the initial particle A after a single reaction A — > B. If it is repeated n 
times in a finite total time T one obtains (using n = J^, where St is the duration of each reaction) 



In the limit of very large n (very small St) we expand the hyperbolic functions in a Taylor series 
and keep terms up to second order in St. We obtain 



P a (\qA,PA>) = ex P(~27t ( pA + qA ^ 2<K 2 + 2^ cosh ^ + 
+ sinh St) 2 - (cosh St + sinh St))) 



(6) 



p n (\q A ,PA>) = exp(-^- t (p A + q A ) 2 (\ + \^ + ^ + 2<*t) - 

St 2 T 
-(! + — + St)) = eM~ m (PA + q A ?St 2 ) = 



(7) 



= exp(- — (pa + qA) 2 St 



T 
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Thus, in the limit n — > oo (St — > 0) we have 




T 



(8) 



n 



Now, although we refer in the former equations to the direction A — > B all our discussion remains 
valid also for the opposite one B — > A. That is, repeating either side of the reaction A <-> .B a large 
number of times n in a finite total time T results, in the limit of very large n, in remaining (with 
probability 1) with the initial particle of the repeated reaction. 

We, now, discuss the reversible reaction A + B <-> C + D in which we have two reacting particles. 
We continue to use the number evolution operator N and take into account that the initial particles 
A and B interact. Thus, representing these particles as the coherent states \qA,PA> and \qB,PB> 
we write, for example, the left hand side direction of the former reversible reaction A + B — > C + D 
as 



where the terms QaQb and PaPb represent, as for the boson particles discussed in (Klauder and 
Sudarshan, 1968), the interaction of the particles A and B, and Na, Nb are the number operators 
for them. Note that, as for the reaction A <-> B (see the discussion after Eqs (4) and (5)), the 
operation of the evolution operator on the coherent state \qA,PA > \qb,Pb >, which is now more 
complicated due to the interaction between A and B, is represented by \qciPc > \QDiPD >■ We 
calculate, now, the probability that the reaction A + B — > C + D results in remaining with the 
initial particles A and B only (we denote this probability by p(\qB,PB> \qa,PA>))- 



Using Eqs (1), (3) and the following coherent states properties (Klauder and Sudarshan, 1968) 
< QiP\Q\liP >= 9) < Q,P\P\l,P >= P (derived by using the operator U from Eq (2) and the 
relation iV|0,0>= 0) we obtain 



exp((iV A + N B + PaPb + QAQB)t)\qA,PA> Wb,Pb>= \qc,Pc> Wd,Pd >, 



(9) 



p(\qb,pb> Wa,PA>) = 



(10) 



=<Qb,Pb\ <qA,PA\exp((N A + N B + PaPb + QAQB)t)\qB,PB> Wa,PA> 



p(\lB,PB> \qa,PA>) = exp((q A qB + PAPB)t) ■ 



• <Qb,Pb\ <QA,PA\qB t ,PB t > \qA t ,PA t >= ex.p((q A qB + PAPB)t) ■ 
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■exp(--(£4 +pa) 2 - ^(qb +pb) 2 - -^{qA t + PA t f - -^{qB t +PB t ?) ■ 

m '" {qA+PA) m {qA t +PA t T , . > S ' r ^° {qA+PA) s {qA t +PA t y 
mjn=0 2 2 ( m !n!)2 s>r=0 2 2 ( s !r!)a 

• <s|r>= exp((g j4 QB + PAPB)t)exp(-^(q A + p A ) 2 - ^(q B + Pb) 2 ~ 

- ^ +pb ( ) ) 2. ( } 

n=0 

" 2^ = exp((gA9B +p J 4Ps)t)exp(--(g j4 - 

r=0 

-^0?B +PB) 2 - ^(<7At +PA t ) 2 ~ ^{qBt+PBtf + \{qA+PA){qA t +PA t ) + 

+ ^(<7B +PB){qB t +PB t )) = exp(^(cosht + smht)((q A + p A ) 2 + (q B + Pb) 2 ) - 
-^(1 + (cosht + sinht) 2 ) • ((q A + Pa) 2 + (qB+PA) 2 ) + {qAqB +PAPbY) 



This is the probability to remain with the original particles A and B after one reaction. Repeating 

T 

St 



it a large number of times n in a finite total time T, where n = £ (5t is the time duration of one 



reaction) one obtains for the probability to remain with A and B. 

T 

P n (\qB,PB> Wa,PA>) = ex.p(—((q A qB + PAPB)5t + 

+ (i(cosh<ft + sinh5t) - i(l + (cosh St + sinh 5t) 2 )) ■ (12) 
i(qA + PA) 2 + (q B + Pb) 2 ))) 

In the limit of very large n we expand the hyperbolic functions in a Taylor series and retain terms 
up to second order in St. Thus, 

T 1 St 2 

P n (\qB,PB> \qA,PA>) = exp(—((q A q B + p A p B )St + (-(1 + — + 

+St)-^(2 + 25t 2 + 25t))((q A +PA) 2 + (q B +PB?))) = (13) 

= exp(T((gA<7B +PaPb) ~ -j((qA +PA) 2 + (qB + Pb) 2 ))) 

Taking the limit of n — ► 00 (St — > 0) we obtain 

lim P n (|g B ,pB> |g J 4,PA>) = exp(T(qAqB + PAPb)) (14) 
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The last probability tends to unity when the c-numbers of either A or B (or both) are zeroes, that 
is, when A or B (or A and B) are in their ground states (in which case they are represented by only 
the first term of the sums in Eqs (1)). Needless to say that all the former discussion remains valid 
also for the opposite direction A + C — > A + B. Thus, we conclude that when either direction of the 
reversible reaction A + B <^ C + D is repeated a large number of times n in a finite total time and 
when at least one of the reacting particles was in the ground state so that its c-numbers are zeroes 
one obtains, in the limit of very large n, a result as if the repeated reaction did not occur at all. 

It can be shown that the general reversible reaction A-y + A 2 + • • • A r <-> B\ + B 2 + ■ ■ ■ + B s , 
where r, s are any two positive natural numbers, also results in a similar outcome if at least one of 
the reacting particles has zero c numbers. We note that the last condition is not necessary when we 
begin with only one reacting particle as we see from the reaction A <-> B. 



3. The probability to find given final configuration different from the initial one 

We now discuss the more general and natural case in which we have an ensemble of particles and 
we calculate the probability to find at the time t a subsystem of this ensemble at some given 
configuration if at the initial time to it was at another prescribed configuration. We assume that the 
corresponding time difference (t — to) is not infinitesimal and that during this time the subsystem 
has undergones a series of reactions. The passage from some reaction at some intermediate time U to 
the neighbouring one at the time (tj + 5t) is governed by the correlation between the corresponding 
resulting configurations of the subsystem at these times. Thus, restricting, for the moment, our 
attention to the case in which a particle in the subsystem that was at the time to i n the state A 
and at the time (to + St) in B we can write the relevant correlation (Klauder and Sudarshan, 1968) 
between these two states as 

T(A,B;to A ,(t + 6t) B ) =<qA tQ ,PA t0 \qB t0+H ,PB t0+6t >, (15) 

where \qA t() iPA t{) > an d \qB to+St ,PB to+St > are the coherent states that represent the particles A and B 
at the times to and to + 5t respectively (see Eqs (3), (4)) and the angular brackets denote an ensemble 
average over all the particles of it. We note that if A = B, r measures (Klauder and Sudarshan, 
1968) the autocorrelation of either the particle A or B, and when A / B r is the crosscorelation 
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(Klauder and Sudarshan, 1968) of the two particles. It can be shown, using Eqs (1) and (15) that 
the following relation 



\t(A, B; to A , (t + St) B )\ 2 = r(A, A; to A , to A )r(B, B; (t + 8t) B , (to + St) B ), (16) 

is valid. That is, the modulus of the crosscorrelation of the particles A and B at the times to and 
(to + St) equals the product of the autocorrelation of the particle A at the time to by that of B at 
the time (to + St). The last relation is interpreted (Klauder and Sudarshan, 1968) as the probability 
density for the occurence of the reaction A — ► B at the time (to + St). That is, given that the system 
was in "state" A at the time to, the probability to find it at the time (to + St) in "state" B is given by 
Eq (16). We can generalize to the joint probability density for the occurence of n different reactions 
between the initial and final times to and t, where each involves two different particles and is of the 
kind Ai — > A,i + \. That is, each reaction is composed of two parts; the first one is that in which a 
particle of the subsystem is observed at the time U to be in state Ai, and the second that in which 
it is observed at the time U + St to be in the state Ai + ±. Thus, the total time interval (t — to) is 
partitioned into 2n subintervals during which the n reactions occur. The total correlation is 

\t(A x ,A 2 ,--- ,A 2n ;t ,t + St,--- ,t)| 2 = T(A 1 ,A 1 ;t ,to) ■ 

■t(A 2 , A 2 ; t + St, to + St)--- r(A 2n , A 2n ; t, t) = (17) 

k=2n— 1 k=n— 1 

= Yl T(A k+1 ,A k+1 ;t + kSt,t + kSt) = J| r(A 2k+1 , A 2k+1 ; t + 

k=0 k=0 

+2kSt, t + 2kSt) ■ r(A 2k+2 , A 2k+2 ;to + (2k + l)5t, t + (2k + l)St) = 

fc=n— 1 

= II |r(A 2fc+1 ,^ 2fc+2 ;t + 2Mt,to + (2A; + l) ( 5t)| 2 

k=0 

The last result was obtained by using Eq (16). By the notation A 2n we mean, as remarked, that 
there are n separate reactions each of which involves two states (and not 2n different particles). 
Now, we show in the previous section, for either direction of the reversible reaction A <-> B, that 
the probability to remain in the initial state A (or B) tends to unity in the limit of a very large 
number of repetitions, in a finite total time, of A — > B (or B — > A) which amounts to performing 
each such reaction in an infinitesimal time 5t. That is, in this limit of vanishing St each factor of 
the last product in Eq (17), which is the probability for the reaction Ai — > Ai + \, tends to unity and 
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with it the joint probability for the occurence of the n reactions. Thus, the specific prescribed path 
of reactions is followed with a probability of unity. 

From the last discussion we may obtain the joint probability density for the case in which some 
of the n intermediate reactions are of the more general kind A\ + A2 + • • • + A r — > B\ + B 2 + ■ ■ ■ + B r , 
where r is an arbitrary natural positive number. That is, at some of the 2n times there may occur, 
in a simultaneous manner, r different reactions at r different places each of the kind A — ► B. 
Thus, we assume that r particles in the subsystem that were at the time to + iSt in the states 
Aj (j = 1, 3, 5, • • • , 2r — 1) were observed at the time (to + (i + l)St) to be in the states Aj+i 
(j + 1 = 2, 4, 6, • • • , 2r). We assume that at each of the other intermediate times there occurs only 
one single reaction Ai — > Ai + \. Thus, there are (n + r — 1) reactions each of them occurs between 
two particles. In this case the corresponding total coherence among all these reactions is 

Ttotal = T(Ai, A 2 , ■ ■ ■ , A i+ i,A i+2 , • • • , A i+2 r, ■ ■■ , -A(2ra+2r- 2) ! *0j *0 + St, ■ ■ ■ 

• • • , to + iSt, t + {i + l)St, • • • , t + iSt, t + (i + l)St, ■■■,*), 

S v ' 

r 

where the underbrace denotes that the r particles observed at the time (to + iSt) as Aj were seen 
to be at the time to + (i + l)St as Aj + \ (j = 1, 3, ■ ■ ■ , 2r — 1). Again the notation A( 2n +2r-2) means 
that we have (n + r—1) reactions each involving, as remarked, two states. Using Eqs (16)-(17) and 
the former equation for Ttotal we find that the joint probability to find at the time t the relevant 
subsystem at the given final configuration after the occurence of these (n + r — 1) reactions is 

\notai\ 2 = t(A ± , A 1 ;t ,t )T(A2, A 2 ;t + St,t + St) ■ ■ ■ r(A i+1 , A i+1 ;t + 
+iSt, t + iSt) • T(A i+2 , A i+2 ; t + (i + l)St, t + (i + l)6t) • • • 

• • • T(A i+2r -l, A i+2r -i; to + iSt, t + iSt) ■ T(A i+2 r, A i+ 2r, to + (« + l)St, t + 

+ (i + l)St) • • • T(A(2n+2r-2),A( 2n+ 2r-2),t0 + ^St, t + 2nSt) = 
k=n— 1 

= II T ( A 2k+i,A 2 k+i;t + 2k5t,t + 2kSt)r(A 2 k+2, A 2 k + 2;t + 
fc=o 

r-l 

+(2/e + l)St, t + (2A; + l)St) ■ JJ r(^ i+i , t + z<5t, t + i<5t) • (18) 

fc=n— 1 

■r(Aj +J - + i, A i+j+1 ;t + (i + l)<5t,t + (i + l)St) = JJ |r(^ 2 fc+i, ^2fc+2; *o + 

fc=0 
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j=r- 1 

+2kSt,t + {2k + l)5t)\ 2 [J |r(A i+j ,^ 

3=1 

The first product of the last result is the same as that of Eq (17) and the second takes account of 
the r — 1 simultaneous reactions at the time (to + (i + l)6t) (the first product involves also one of 
the r simultaneous reactions at the time (to + (i + 1)<^))- Each of the reactions in both products is 
of the kind A — > B which was shown in the former section (see also the discussion after Eq (17)) 
to have a unity probability for remaining with the initial particle A in the limit in which the time 
duration of the reaction becomes infinitesimal. That is, in this limit in which the time alloted for 
each reaction A — > B becomes very small each factor of each product of Eq (18), and with it the 
whole expression, tends to unity. If any particle A of the subsystem does not react with any other 
particle at some of the 2n intermediate times then we may denote its no-reaction at these times as 
A — > A and the probability for it to remain in the initial state (which is the same as the final one) 
is obviously unity. 

Thus, we see that the probability to find at the time t the given ensemble of particles following 
a given path of reactions (from a large number of possible paths) tends to unity if the relevant 
reactions are performed in a dense manner. That is, each occuring in an infinitesimal time 5t. 

4. Billiard simulation of the reversible reaction A + B <-> A + C 

We, now, show that the results described in the former sections have also strong numerical support. 
This is demonstrated for the reversible reaction A + B <-> A + C which is simulated here by using the 
two-dimensional circular billiard (Bar, 2001) which is composed of two concentric circles. We assume 
that initialy we have a large ensemble of identical point particles each of them is the component A 
of the reaction A + B <-> A + C . All of these particles are entered, one at a time, into the billiard in 
which they are elastically reflected by the two concentric circles. That is, the angles before and after 
each reflection are equal. We assume that on the outer circle there is a narrow hole through which 
the particles A leave the billiard. Once a particle is ejected out a new one is entered and reflected 
inside the billiard untill it leaves and so for all the particles of the ensemble. There are two different 
possible kinds of motion for each point particle A before leaving the billiard; either it is reflected 
between the two cocentric circles or, when the angle of reflection is large, reflected by the outer 
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circle only without touching the inner one. Now, since both motions are elastic each particle A, once 
it begins its reflections in either kind of motion, continues to move only in this kind untill it leaves 
through the narrow opening. The component B of the reaction denotes the outer larger circle, and 
the component C denotes both circles. That is, the left hand side A + B of the reaction signifies 
that the point particle A moves inside the billiard and is reflected by the outer circle only, whereas, 
the right hand side A + C denotes the second kind of motion in which the point balls A are reflected 
between the two circles. We call these two kinds of paths "states" (Bar, 2001), so that the path that 
touches both circles is "state" 1 and the one that thouches the outer circle only is "state" 2. This 
billiard model was studied in (Bar, 2001) as an example of a classical system that behaves the same 
way quantum systems do when exposed to a large number of repetitions, in a finite total time, of 
the same experiment ((Misra and Sudarshan, 1977), (Giulini et al, 1996), (Pascazio and Namiki, 
1994), (Cook, 1988), (Harris and Stodolsky, 1981), (Bixon, 1982), (Itano et al, 1990), (Kofman and 
Kurizki, 1996), (Simonius, 1978), (Aharonov and Vardi, 1980), (Facchi et al, 1999)). The concentric 
billiard is, schematically, represented in Figure 1. 




Figure 1. A schematic representation of the concentric circular billiard that simulates the reversible reaction 
A + B <-> C + D. 

Now, since in such a system we can not follow the path of each particle and can not differentiate 
between the two kinds of motion we have to consider, as done for the nuclear and radioactive 
processes (Bauer and Bertsch, 1990), the activities of these particles in either path. That is, the rate 
at which the entire ensemble of particles, being at either state, leaves the billiard. We assume for the 
activity discussed here, as is assumed (Bauer and Bertsch, 1990) for the nuclear and radioactive's 
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activities, that each particle A enjoys arbitrary initial conditions, so in the following numerical 
simulations we assume that it may begin its journey inside the billiard at either "state" which is 
determined randomly using a random number generator. 

We want to show, numerically, that if either side of this reversible reaction A + B ^ A + C is 
repeated a large number of times N in a finite total time T, then, in the limit of very large N , 
the activity obtained is the same as the "natural activity" (Bar, 2001) that results when no such 
repetitions are performed (Bar, 2001). For that matter, we take into account that the reversible 
reactions that occur in nature have either equal or different rates for the two directions of the 
reactions and that the total activity of the ensemble of particles depends critically upon these rates 
(Bauer and Bertsch, 1990). If, for example, we consider the equal rate case then we have to discuss 
the rate of evacuation of the billiard when each particle is allowed, after a prefixed number of 
reflections in either state, to pass, if it is still in the billiard, to the other one. This activity is shown 
by the solid curve in Figure 2 in which the ordinate axis denotes the number of particles A that 
leave the billiard in prescribed time intervals binned in units of 60 (Bauer and Bertsch, 1990). We 
assume (Bauer and Bertsch, 1990) that each point particle A in either state moves with the same 
speed of 3, and the hole through which they leave has a width of 0.15. We denote the outer and 
inner radii of the billiard by r\ and r2 respectively, and assign them the values of r\ = 6 and r<i = 3. 
The initial number of the particles A was 10 6 , and each one of them passes from one "state" to the 
other, if it did not leave the billiard through the hole, after every 1100 consecutive reflections. We 
note that this rate of one passage for every 1100 reflections is typical and common for these kinds 
of billiard simulations (Bauer and Bertsch, 1990), (Gutkin, 1986), (Hobson, 1975). 

The natural activity is obtained, as remarked, when the entire ensemble of 10 6 particles A enter, 
one at a time, to the billiard at the same "state" and remain all the time in this "state" without 
passing to the other until they leave the billiard. The dashed curve in Figure 2 shows this natural 
activity when all the particles A are in "state" 2 in which they are reflected only between points of 
the outer circle untill they leave the billiard. The dash-dot curve shows the activity when all the 
particles A are in "state" 1 in which they are reflected only between the two circles. It has been 
found that for the values assigned here to the radii of the outer and inner circles (6 and 3) the 
activity of "state" 2 shown by the dashed curve is the maximum avilable and that of state 1 shown 
by the dashdot curve is the minimum. The large difference between the two activities has its source 
in the range of the allowed angles of reflections which is much larger in state 2 than in state 1. 
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This is because the minimum trajectory between two neighbouring reflections in state 2, where the 
particles A are reflected between points of the outer circle only, may be infinitesimal compared to 
the corresponding trajectory in state 1 which is (we denote the trajectories between neighbouring 
reflections in states 1 and 2 by d\ and d 2 respectively) d\ min = r\ — r 2 . For the values assigned here 
to the radii r± and r 2 of the two concentric circles (ri = 6 and r 2 = 3) d\ min = 3. We note that the 
maximum trajectory between two neighbouring reflections in state 2 is equal to the corresponding 
one in state 1, that is 

di max = d 2max = \jr\ + rl 

Thus, the particles in A have many more possibilities to be reflected to the hole and leave the billiard 
in state 2 than in state 1 and, accordingly, their activity is much larger. The solid curve in Figure 2 
is, as remarked, the activity obtained when the particles A are transferred between the two states 
at the rate of one passage for every 1100 reflections and so, as expected, its activity is between the 
two other activities shown in Figure 2. 




Figure 2. The dashed curve shows the activity obtained when all the 10 particles of the ensemble are allowed to be 
only in "state" 2, in which they are reflected between points on the outer circle only. The dashdot curve is the activity 
when all these particles are allowed to move only in "state" 1 (between the two circles). For the values assigned here 
to the outer and inner circles (6 and 3) the dashed (dashdot) curve is the maximum (minimum) activity. The solid 
curve shows the activity obtained when the particles in either state pass to the other after every 1100 reflections. The 
x axis specifies time binned in units of 60. 

We numerically interfere with the rate of the systematic passage of the point particles A between 
the two states such that this rate is accelerated. We refer in the following not only to a specific 
passage between the two "states" but also gives the results in parentheses for the opposite passage. 
Thus, it is found that the activity of the entire ensemble is directly (inversely) proportional to the 
rate of the passage from state 1 (2) to state 2 (1) when the opposite passage from state 2 (1) to 
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state 1 (2) remains at the rate of one for every 1100 reflections. Thus, we have found that when the 
particles in state 1 (2) are transferred to state 2 (1) at the maximum rate of one passage for each 
single reflection and the particles in state 2 (1) are passed to the state 1 (2) at the rate of one for 
every 1100 reflections then the activity of the particles A is maximal (minimal). But as we have 
remarked the maximal (minimal) activity is obtained only when each particle of the entire ensemble 
is always in state 2 (1). In other words, as we have remarked, a very large number of repetitions 
of the left (right) direction A + B — > A + C (A + C — > A + B) of the reaction where the right 
(left) direction A + C — > A + B (A + B — > A + C) occurs every 1100 reflections, yields a result 
as if the densely repeated reaction never happened and the activity obtained is the natural one in 
which no external repetition is present. The dashed curve in Figure 3, which is the same as the 
dashed one of Figure 2, shows the activity obtained when all the 10 6 particles A of the ensemble are 
allowed to move only in state 2 until they leave the billiard. The solid curve is the activity obtained 
when the reaction A + C A + B is repeated after each single reflection and the opposite one 
A + B — > A + C after every 1100 reflections. It is seen that these curves are similar to each other. 
That is, we realize, in accordance with the analytical results of the former sections, that a large 
number of repetitions of the reaction yields a result that characterizes the activity obtained in the 
absence of such repetitions. 
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Figure 3. The dashed curve, which is the same as the dashed curve of Figure 2 (they look slightly different since the 
abcissa axes of these figures are different), shows the activity obtained when all the 10 6 particles A of the ensemble 
are numerically constrained to be only in "state" 2 until they evacuate the billiard. State 1 is not allowed for them. 
The solid curve is the activity obtained when each particle in state 1 is passed to "state" 2 after each single reflection, 
whereas those in "state" 2 pass to the opposite one only after every 1100 reflections. As for figure 2 the abcissa axis 
denotes time binned in units of 60. Note the similarity between the two curves. 

This is seen, in a much more clear way, in Figure 4 for the other direction A + B^A + Coi the 
reaction. The apparent single graph of the figure is actually composed of two curves; one solid and 
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the other dashed. The solid curve shows the activity obtained when the reaction A + B ^ A + C is 
repeated after each single reflection and the opposite one A+C — ► A+B after every 1100 reflections. 
The dashed curve, which is identical to the dash-dot one from Figure 2, is the activity obtained 
when all the particles A of the ensemble are constrained to move only in state 1 until they leave 
the billiard. Note that the two curves are almost the same except for the longer tail of the dashed 
curve. 
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Figure 4- The apparently one curve shown in the figure is actually two curves one dashed and the other solid. The 
dashed curve, which is the same as the dashdot curve of Figure 2, shows the activity obtained when all the 10 6 
particles A of the ensemble are numerically constrained to be only in "state" 1 until they evacuate the billiard. State 
2 is not allowed for them. The solid curve is the activity obtained when each particle in state 2 is passed to "state" 
1 after each single reflection, whereas those in "state" 1 pass to the opposite one only after every 1100 reflections. As 
for figure 2 the abcissa axis denotes time binned in units of 60. Note that the two curves are almost identical (the 
dashed curve has a longer tail (for large t) than the solid one). 

From both Figures 3 and 4 we realize that the large number of repetitions of either direction of 
the reversible reaction A + B <-> A + C has the effect as if it has not been performed at all and the 
actual activity obtained is that of the natural one that does not involve any repetitions. 

We note that as the analytical results are obtained in the limit of the largest number (actually 
infinite) of repetitions so the similar numerical results are obtained in the limit of the largest number 
of repetitions of the reaction. That is, of numerically repeating it after each single reflection. In other 
words, a mere high rate (which is not the maximal) of one side of the reaction is not enough to 
produce the results shown in Figures 3-4. This is clearly shown in Figure 5 the solid curve of which 
shows the activity obtained when each particle in "state" 1 is passed to "state" 2 after every two 
neighbouring reflections whereas those of "state" 2 are passed (one at a time) after every 1100 
reflections. Note that the solid curve in Figure 3 shows the activity obtained when the particles in 
"state" 1 are passed to "state" 2 after each reflection and those of 2 passed to 1 after every 1100 
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reflections. That is, although the two high rates represented by the two solid curves in Figures 3 
and 5 are almost the same nevertheless the resulting activities, contrary to what one may expect, 
are very different. That is, that of Figure 3 is much higher than that of Figure 5 as may be seen 
from the solid curve that begins at t = 60 (note that our abcissa axis is binned in units of 60) from 
the high value of 5.65 • 10 5 and ends at t = 360. The corresponding solid curve of Figure 5 begins 
at t = 60 at the much smaller value of 4.45 • 10 5 and ends at the later time of t = 420. That is, 
by only increasing the rate of repeating the same reaction from one for every two reflections to one 
for each reflection results in an additional 120000 particles that leave the billiard already at the 
first binned time unit. The two dashed curves of Figures 3 and 5 are identical and denote the same 
activity obtained when all the 10 6 particles A of the ensemble are numerically constrained to be 
only in "state" 2 until they evacuate the billiard. 



5 - 1 
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Figure 5. The dashed curve, which is the same as the dashed curve of Figure 3, shows the activity obtained when all 
the 10 6 particles A of the ensemble are moving inside the billiard only in "state" 2 until they are evacuated outside. 
State 1 is not allowed for them. The solid curve is the activity obtained when each particle in state 1 is passed 
to "state" 2 after every two reflection, whereas those in "state" 2 pass to the opposite one only after every 1100 
reflections. Note that although the two solid curves of Figures 3 and 5 differ slightly in their reaction rates (of one for 
each reflection in Figure 3 and of one for each two reflections in Figure 5) the activities are very different (see text). 
As for all the former figures the abcissa axis denotes time binned in units of 60. 

Thus, as remarked, the important factor that causes a result of maximum activity is the highest 
possible rate and not merely a large ratio between the higher and slower frequencies. This is in 
accord with the analytical results of Sections 2 and 3 which are obtained for the largest rate (actually 
infinite) of repeating any direction of the general reversible reaction A\ + Ai + • • • A r <-> B\ + B2 + 
• • - +B S , where r and s are any two arbitrary natural positive numbers. In this limit we find, as found 
from the billiard simulation, that one remains, with a unity probability, with the initial reacting 
particles as if the repeated reaction did not occur at all. 
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All the former simulations were done when the outer and inner circles radii were 6 and 3 
respectively. We note that we obtain similar numerical results for all other assigned values of r\ 
and T2 up to the extreme limits of r\ » and T\ r2 provided we always have r\ > T2- 

We may explain these results along the same line of explanation used to interpret the similar re- 
sults obtained in the quantum regime. It have been established theoretically ((Misra and Sudarshan, 
1977), (Simonius, 1978), (Aharonov and Vardi, 1980), (Pascazio and Namiki, 1994), (Facchi et al, 
1999)) and experimentally ((Itano et al, 1990), (Kofman and Kurizki, 1996), (Giulini et al, 1996)) 
that taking some quantum system which may reduce by experimenting on it to any of its relevant 
eigenstates and repeat this experiment a large number N of times in a finite total time T results, 
in the limit of very large N (keeping T constant), in preserving the initial state of the system. 
This phenomenon is the Zeno effect ((Misra and Sudarshan, 1977), (Simonius, 1978), (Aharonov 
and Vardi, 1980), (Pascazio and Namiki, 1994), (Facchi et al, 1999)). The similar results obtained 
theoretically in Section 2 suggest that this effect may be effective also in the classical reactions. That 
is, repeating these reactions a large number of times, in a finite total time, may results in remaining 
with the initial reacting particles (initial "state"). Moreover, it have been shown ((Aharonov and 
Vardi, 1980), (Facchi et al, 1999)) that if any quantum system evolves in a number of different 
possible paths of states, each of which begins at the same given initial state and end at another 
given common final one, then it is possible to realize any such path by making dense measurement 
(in a finite total time) along it. That is, by performing in a dense manner the set of experiments 
that reduce the system to the specific states that constitute the relevant path. Similar results were 
obtained in Section 3, in which we show that the joint probability density for the occurence of n 
given different reactions between the initial and final times to and t tends to unity in the limit of 
densely performing these reactions (when n — > oo). 

The same results were obtained also in the billiard simulations from which we realize that 
repeating a large number of times any direction of the reversible reaction A+B <-> A+C have, in the 
limit of numerically repeating it after each single reflection, the effect as if it has never occured. This 
is what one obtains in the Zeno effect ((Misra and Sudarshan, 1977), (Simonius, 1978), (Aharonov 
and Vardi, 1980), (Pascazio and Namiki, 1994), (Facchi et al, 1999), (Itano et al, 1990)) where the 
system is preserved in the initial state in spite of the very large number of measurements. 
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Concluding Remarks 



We have discussed classical reactions using quantum theory methods in which particles are rep- 
resented by coherent state functions and the product of these states with their conjugates is 
interpreted as probability (Masao, 1976). It is shown that if either direction of the general reaction 
A-y + A 2 + • • • A r <-> B\ + B 2 + • • ■ + B s , where r and s are any two natural positive numbers, is 
repeated a large number N of times in a finite total time then in the limit of a very large N one 
remains, with a unity probability, with the initial particles only. In this context we differentiate 
between the case in which there were more than one initial reacting particle and the case in which 
there is only one such particle. In the first case the mentioned unity probability is obtained if at 
least one of the initial reacting particles have zero values for its c-numbers q and p that denote its 
coherent state representation, whereas, in the second case this condition is not needed. Moreover, 
it has been shown that any prescribed evolution (from a large number of possible ones) through a 
sequence of specific reactions may be realized with a unity probability by densely performing these 
reactions. This effect that results from increasing the rate of the reaction have been demonstrated 
also through numerical simulation for the reversible reaction A+B <-> A+C. We use for that purpose 
the two-dimensional concentric billiard in which the two modes of possible reflections represent the 
two sides of the reaction. It has been shown that by repeating either side of the reaction a large 
number N of times (whereas the other side was repeated with a much lower rate) we obtain, in the 
limit of the largest rate of repetition, a result as if the fastly repeated reaction were not performed 
at all. 

The obtained results conform to the Zeno effect which is considered in (Simonius, 1978) to hold 
also in classical and macroscopic phenomena. In this effect ((Misra and Sudarshan, 1977), (Simonius, 
1978), (Pascazio and Namiki, 1994), (Itano et al, 1990)) the very large number of repeating the same 
experiment, in a finite total time, results in preserving the system in the state it was before initiating 
these repetitions. Moreover, it has been shown in ((Aharonov and Vardi, 1980), (Facchi et al, 1999)) 
that this effect can be generalized to a whole path of states in which the final state is different 
from the initial one. That is, the mechanism of dense measurement causes the "realization" of this 
specific path from a very large number of different ones. This result was obtained by calculating the 
joint probability of n different reactions, each of which occurs at its specific place and time, where 
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we see that densely performing these reactions causes the joint probability for the occurence of the 
n reactions (where n — ► oo) to tend to the unity value. 
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